X-point collapse and saturation in the nonlinear tearing mode reconnection 
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We study the nonlinear evolution of the resistive tearing mode in slab geometry in two dimensions. 
We show that, in the strongly driven regime (large A'), a collapse of the X-point occurs once the 
island width exceeds a certain critical value ~ 1/A'. A current sheet is formed and the reconnection 
is exponential in time with a growth rate oc rj^^^, where 77 is the resistivity. If the aspect ratio of 
the current sheet is sufficiently large, the sheet can itself become tearing-mode unstable, giving rise 
to secondary islands, which then coalesce with the original island. The saturated state depends on 
the value of A'. For small A', the saturation amplitude is oc A' and quantitatively agrees with the 
theoretical prediction. If A' is large enough for the A-point collapse to have occured, the saturation 
amplitude increases noticeably and becomes independent of A'. 

PACS numbers: 52.35. Vd, 52.65.Kj, 52.35.Py 



Magnetic reconnection is the breaking and rejoining of 
magnetic field lines in a plasma. Solar flares are believed 
to be a manifestation of this phenomenon . Other clas- 
sical examples are reconnection between the solar and 
the Earth's magnetic field in the magnetopause and the 
magnetotail and the sawtooth instability in tokamaks 
0. In some cases, most notably the sawtooth, recon- 
nection takes place in a plane perpendicular to a strong 
magnetic field, in which case it occurs via the tearing- 
mode instability. Linear theory shows that an MHD 
equilibrium is tearing-mode unstable if the instability pa- 
rameter A' > 0. Analytical and numerical studies of the 
tearing mode have been mostly restricted to low values of 
A'. However, it has been shown that kinetic effects can 
change the instability threshold to A' > A'^^,^^. > 1 0, Q 
and there is, indeed, experimental evidence for A' 3> 1 in 
the sawtooth 0. The evolution of large-A' (i.e., strongly 
driven) tearing modes, even in the simplest physical mod- 
els, remains poorly understood. To address this problem, 
we investigate the evolution of the tearing mode in the 
broadest ranges of A' and the resistivity rj achieved to 
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FIG. 1: Effective growth rate at the A-point 7cff = dln^/df 
vs. time for a strongly driven (large A') tearing mode. 



date. We find that, for sufficiently large A' and suffi- 
ciently small 7], the tearing modegoes through five stages 
(Fig. nj : (I) linear instability |J| , (II) algebraic growth 
(Rutherford stage) . (Ill) AT-point collapse followed by 
current-sheet reconnection (Sweet-Parker 0,01 stage), 
(IV) tearing instability of the current sheet resulting in 
generation of secondary islands, and (V) saturation. The 
traditional theory of the tearing mode, valid for small A', 
comprises just Stages I, II and V. The fact that, at large 
A', Stage II is succeeded by Stages III, IV, and a modified 
Stage V is new in the tearing- mode context. 
We solve the Reduced MHD equations 



dtoj + v± ■ Vuj = B± ■ Vj|| , 



(1) 
(2) 



in a two-dimensional periodic box x Ly using a pseudo- 
spectral code at resolutions up to 3072 x 4096. The 
total magnetic field is -B = BzBz + B±, the in-plane 
magnetic field is B± = x 'Vip, the in-plane veloc- 
ity is v± = Ez X V0, and w = • (V x v±) = V'^ij), 
j\\ — • (V X B) = V^V'- We impose the equilibrium 
configuration = t/'o/ cosh^(a;) and = 0. By set- 
ting -00 = 3'\/3/4, we scale the units of field strength in 
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1. All lengths are scaled so that = 27r. 
Time is, therefore, scaled by the in-plane Alfvcn time 
Lx/^irBy^laax- To the equilibrium, we add an initial per- 
turbation -0^^-' = ipi{x) cos(fcy), where k = Lx/ Ly. Given 
a perturbation in this form, the island width W and the 
reconnected flux ^[t) = —ip{t, 0, 0) -I- ipo are related by 



W = 4^*(t)/V^,'(0). (3) 
For our equilibrium, the instability parameter is [l^ 
V'i(+0) - ^i(-O) _ 2(5 - fc2)(3 + fc2) 
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FIG. 2: Contours of '0 at the beginning and end of Stage III 
in Fig.Q The boundaries of these plots are noi the boundaries 
of the computational box. 



The equilibrium is tearing-unstable if A' > A: < -v/S- 
A' is varied by changing k, i.e., Ly. 

We now describe the evolution of the tearing mode 
stage by stage. During Stages I-II, reconnection occurs 
via an X-point configuration. In Stage I, it is a linear in- 
stability with the island width W growing exponentially 
in time Once W exceeds the resistive scale, this stage 
gives way to the Rutherford stage (Stage II), during 
which the growth is algebraic in time: dW/dt ~ 77A'. 
Omitting further discussion of these stages, which have 
been studied before proceed to 

Stage III: X-Point Collapse and Sweet-Parker Recon- 
nection. In simulations with large A', the X-point con- 
figuration eventually collapses and a current sheet is 
formed (Fig. accompanied by a dramatic speed up 
of the island growth (Fig. Several previous numeri- 
cal studies in various reconnection contexts have also re- 
ported a nonlinear speed-up [l^ and a tendency for the 
current-sheet formation Waelbroeck |0| predicted 
that when the island width W > Wc ~ 1/A', no equi- 
librium X-point configuration exists and a current sheet 
must form. By varying rj and A' in our simulations, we 
have tried to verify this prediction. We define Wc as the 
island width at which d'jcff/dt = after the Rutherford- 
like algebraic stage (e.g., at t « 315 in Fig. In Fig. O 
we plot A'Wc vs. 77 for two different values of A'. The 
dependence of A'Wc on 77 appears to be linear and ex- 
trapolates in the limit of 77 to A'Wc — 8.2 for both 
values of A' used. Thus, the transition criterion is 
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FIG. 3: The critical island width for collapse vs. 77 at fixed 
A' — 17.3, 30.1. Dashed lines are linear fits. 



itatively in the way described in |^ ^3 but with the 
outflow velocity 7;out and the current sheet length Leg 
growing with time. Since the reconnected flux '5 changes 
at the SP rate oc 7/^/^, we can assume that the evolu- 
tion is quasistatic, so that the system passes through a 
sequence of ideal equilibria, in each of which L^g and 
the configuration outside (but not inside) the current 
sheet arc fully determined by the instantaneous value of 
Let us assume that in these equilibria, the vicin- 
ity of the current sheet is described by t he Syrovatskii 
solution with a unidirectional current |l7l |. In this solu- 
tion, the magnetic field immediately outside the sheet is 
Bin = Bo{^)Lcs{^)/Lx, where 4' = ^'(i) is the recon- 
nected flux and Bq is the fleld away from the sheet. Then 
the reconnected flux grows according to (cf . [isl ) 



— ~ Win Sin 

dt 
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(6) 



where we have used the SP expression for the inflow ve- 
locity, Win ^ (?7i'out/-^cs)^^^j and taken the outflow veloc- 
ity to be Alfvenic, Vout ~ ^in- Eq. ijHJ) implies that the 
growth of must speed up compared to 5* oc (?/A't)^ 
in the Rutherford [g stage (Stage II). When 5* is close 
to its value ^'c at the beginning of the collapse, we may 
approximate -Bo(^) ~ ^0(^*0) = const. This implies 



A'Wc ~ 8.2 + /(A')77, 



(5) 



where the slope function /(A') remains undetermined 
but is seen in Fig. O to increase with A'. 

Fig. Ela) shows that, in this stage, the reconnected 
flux (measured at x ~ y ^ 0) grows exponentially in 
time: ln(^ — V&c) = 7sp(i — tc), where tc is the time at 
which the c ollap se begins, ^'c = ^(^c), and 7sp is the 
growth rate [26j. Varying A', we have ascertained that 
7SP is independent of A'. Its dependence on rj is plotted 
in Fig.^b). The scafing 7sp oc 77^/^ appears to hold. 

We think that what we observe is an exponential-in- 
time Sweet-Parker (SP) reconnection that proceeds qual- 



A'= 17.3 

7=6.5x10-* 
. 7 = 5.0x10-* 
. (, = 2.8x10-' 
- (,= 1.0x10-* 



0.08 
0.07 
0.06 
0.05 
0.04 
: 0.03 



100 



150 



(a) 



(b) 



FIG. 4: (a) Growth of the reconnected fiux "l/ during the SP 
stage for fixed A' = 17.3 and various values of 77. (b) Slopes 
of these lines vs. 77 during the exponential growth. 



3 



1.0t 
0.8- 



0.0 0.2 0.4 0.6 



(a) 



(b) 



FIG. 5: The exponential stage (370 < t < 450) of the run 
of Fig. (a) maximum outflow velocity Vout vs. the current- 
sheet length Lcs; (b) Lcs vs. (^— \l'c)/Bo, where Bq is defined 
as the maximum value of By along the x axis. These curves 
do not extrapolate to the origin because the full-width-half- 
maximum definition used for Lcs correctly reflects the growth 
of the current-sheet length but not its true length (thus, it 
formally gives Lcs > for the X-point reconnection) . 



Wout/-^cs ~ Bin/ Lcs ~ Bo/Lx = const, a conclusion 
confirmed by Fig.[S^a). Since Lcs(^c) = 0, Lcs should 
be a growing function of — ^fc. Indeed, Eq. (j^J) is 
consistent with the numerically observed exponential SP 
reconnection if Lcs ^ - ^c)/Bo [cf. Fig.lHJb)]. 

The elongation of the current sheet ceases when Lcs 
reaches a significant fraction of the box size. Reconnec- 
tion can still proceed in a SP fashion, but the growth of 
the reconnected fiux slows down (see Fig. Indeed, in 
the right-hand side of Eq. (j^J), Lcs no longer increases 
with 5" and i?o(^) starts to decrease as the initial recon- 
nectable fiux is used up. In Fig. we show the current- 
sheet length Lcs and width 6cs measured using a full- 
width-half-maximum estimate at the time when the max- 
imum Lcs is reached. We see that, for fixed A', Lcs is 
roughly independent of r], while Scs ^ V^^'^, in agreement 
with the SP prediction. On the other hand, for fixed rj, 
both Lcs and Sqs grow linearly with A' (cf. except 
for the largest data point, A' = 52.7 

Stage IV: Secondary Island Generation. When the 
aspect ratio of the current sheet A = Lcs/^cs 50, 
the sheet itself becomes unstable to tearing modes and 
generates secondary islands. We expect that this crit- 
ical value is independent of either A' or rj, but due to 





A' = 


30.1 


■ 


A' ^ 


17.3 


• ^cs 


A' = 


30,1 


'■^ 5cs 


A' - 


17.3 



0.08 
0.06 
04 " 



10 

(a) 



10 20 30 40 50 



(b) 




FIG. 6: The current sheet length Lcs and width Scs (a) vs. rj 
and (b) vs. A'. 



FIG. 7: Contours of showing the current sheet instability 
(a-c) and the subsequent nonlinear evolution of the secondary 
island (d-f) for a run with A' = 40.6, r? = 2.8 x 10"''. 



resolution constraints, we do not yet have a numerical 
confirmation of this conjecture [2^ . 

A detailed view of the instability is given in Fig. [7| 
As suggested in [i^, a secondary island first appears 
as a long-wavelength linear perturbation to the current 
sheet, with two X points forming at the ends of the sheet 
[Fig.EJb)]. The reconnection proceeds in a manner anal- 
ogous to Stages I-III discussed above: as the secondary 
island grows, the two secondary X points collapse, giving 
rise to two current sheets, while the island is circularized 
[Fig.EJc)]. The primary and the secondary islands exert 
attracting forces on each other. When the secondary is- 
land is sufficiently large, this attraction causes it to split 
into two parts, which then coalesce with the primary is- 
land [Fig. d-f)]. Note that the sphtting of the secondary 
island into two is a result of the exact symmetry of our 
configuration about the x axis. Even a slight breaking of 
this symmetry would cause the entire secondary island 
to move either upwards or downwards to coalesce with 
the primary |29j . 

We note that, given small enough r], the secondary cur- 
rent sheets should be unstable to generation of tertiary 
islands etc. Also, if the initial flux is not yet exhausted 
after the secondary island has coalesced with the primary, 
the primary current sheet can be regenerated via a second 
collapse (Fig.m Stage IVb). Given sufficiently large A', 
the cycle of current-sheet formation — secondary-island 
generation — coalescence may be repeated several times 
before saturation is reached. 

Stage V: Saturation. The saturated island width in 
the limit of small A' has recently been calculated by Es- 
cande & Ottaviani and Militello & Porcelli a 
theory henceforth referred to, using a liberal permutation 
of the first letters of the authors' surnames, as POEM. 
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FIG. 8: (a) Saturated amplitude ^aat vs. A' for different 
values of r/. The theoretical curves by POEM [Eq. l|7|] and 
White et al. are also shown. The island width formula Q 
has been used to convert V14at calculated by these authors into 
v&aat. (b) *sat VS. Tj foT A' = 8.2, 17.3. In both plots, hollow 
points are the cases where Wsat exceeded the box size. 



They found 



2.44a2A', 



a' = -V'(,'(0)/C'(0). (7) 



For our equilibrium, = 0.125. Fig.|Slshows the depen- 
dence of the numericaUy obtained saturated flux on A' 
and rj compared to the quantitative predictions of POEM 
and of the earlier theory of White et al. [i^. We plot 
^fsat instead of M4at because, for the largest A' values, 
the island width exceeds the box size (in which case 
the saturation is likely to be strongly dependent on the 
equilibrium configuration). For A' < 5, there is excellent 
agreement with POEM [Eq. iQ], but not with White 
et al. [23. The occurence of the X-point collapse, i.e., 
whether the saturation is achieved via current-sheet or 
X-point reconnection, changes the saturated state: Fig.|Sl 
shows a jump in ^'sat at values of A' and ry for which the 
X-point collapse took place in Stage III. For sufficiently 
small 77, the saturated amplitude does not depend on rj 
[Fig. IHIb)]. Also, 'I'sat appears to reach a plateau for 
large A' [Fig. |HIa)] , so that Wsat ~ system size. 

Note that the collapse can occur only if the saturated 
island width is larger than Waelbroek's critical value 
(Fig. EJ, Wsat > Wc. Using Eq. © and Eq. 0, this 
gives A' > 5.2 in the limit of 77 ^ 0. 

In this Letter, we have shown that, for a simple resis- 
tive model of the tearing mode with large A', X-point 
reconnection gives way to much faster current-sheet re- 
connection. Furthermore, the current sheet cannot ex- 
ceed a certain critical aspect ratio, above which it frag- 
ments into secondary islands and current sheets. We 
believe the rather complex behavior we have identified 
to be a generic feature of strongly driven reconnection. 
However, a caveat is in order. While the large-A' con- 
figurations are often encountered in laboratory reconnec- 
tion, understanding the physics responsible for setting 
up these configurations remains a theoretical challenge. 
This unknown physics, along with a number of kinetic 
effects known to be im por tant in various laboratory and 
astrophysical contexts |2J|, must, strictly speaking, be a 



part of any quantitative description of the tearing-mode 
reconnection in real plasmas. 
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